\(\int \frac {(A+B x) (d+e x)^{3/2}}{(b x+c x^2)^3} \, dx\) [1250]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (verified)
   Maple [A] (verified)
   Fricas [B] (verification not implemented)
   Sympy [F(-1)]
   Maxima [F(-2)]
   Giac [A] (verification not implemented)
   Mupad [B] (verification not implemented)

Optimal result

Integrand size = 26, antiderivative size = 346 \[ \int \frac {(A+B x) (d+e x)^{3/2}}{\left (b x+c x^2\right )^3} \, dx=-\frac {\sqrt {d+e x} \left (A b c d+\left (2 A c^2 d+b^2 B e-b c (B d+A e)\right ) x\right )}{2 b^2 c \left (b x+c x^2\right )^2}-\frac {\sqrt {d+e x} \left (b (c d-b e) \left (6 b B c d-12 A c^2 d-2 b^2 B e+7 A b c e\right )-3 c (c d-b e) \left (8 A c^2 d+b^2 B e-4 b c (B d+A e)\right ) x\right )}{4 b^4 c (c d-b e) \left (b x+c x^2\right )}-\frac {3 \left (16 A c^2 d^2+b^2 e (4 B d+A e)-4 b c d (2 B d+3 A e)\right ) \text {arctanh}\left (\frac {\sqrt {d+e x}}{\sqrt {d}}\right )}{4 b^5 \sqrt {d}}+\frac {3 \left (16 A c^3 d^2-b^3 B e^2-4 b c^2 d (2 B d+5 A e)+b^2 c e (8 B d+5 A e)\right ) \text {arctanh}\left (\frac {\sqrt {c} \sqrt {d+e x}}{\sqrt {c d-b e}}\right )}{4 b^5 \sqrt {c} \sqrt {c d-b e}} \]

[Out]

-3/4*(16*A*c^2*d^2+b^2*e*(A*e+4*B*d)-4*b*c*d*(3*A*e+2*B*d))*arctanh((e*x+d)^(1/2)/d^(1/2))/b^5/d^(1/2)+3/4*(16
*A*c^3*d^2-b^3*B*e^2-4*b*c^2*d*(5*A*e+2*B*d)+b^2*c*e*(5*A*e+8*B*d))*arctanh(c^(1/2)*(e*x+d)^(1/2)/(-b*e+c*d)^(
1/2))/b^5/c^(1/2)/(-b*e+c*d)^(1/2)-1/2*(A*b*c*d+(2*A*c^2*d+b^2*B*e-b*c*(A*e+B*d))*x)*(e*x+d)^(1/2)/b^2/c/(c*x^
2+b*x)^2-1/4*(b*(-b*e+c*d)*(7*A*b*c*e-12*A*c^2*d-2*B*b^2*e+6*B*b*c*d)-3*c*(-b*e+c*d)*(8*A*c^2*d+b^2*B*e-4*b*c*
(A*e+B*d))*x)*(e*x+d)^(1/2)/b^4/c/(-b*e+c*d)/(c*x^2+b*x)

Rubi [A] (verified)

Time = 0.58 (sec) , antiderivative size = 346, normalized size of antiderivative = 1.00, number of steps used = 6, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.192, Rules used = {832, 836, 840, 1180, 214} \[ \int \frac {(A+B x) (d+e x)^{3/2}}{\left (b x+c x^2\right )^3} \, dx=-\frac {3 \text {arctanh}\left (\frac {\sqrt {d+e x}}{\sqrt {d}}\right ) \left (b^2 e (A e+4 B d)-4 b c d (3 A e+2 B d)+16 A c^2 d^2\right )}{4 b^5 \sqrt {d}}+\frac {3 \left (b^2 c e (5 A e+8 B d)-4 b c^2 d (5 A e+2 B d)+16 A c^3 d^2+b^3 (-B) e^2\right ) \text {arctanh}\left (\frac {\sqrt {c} \sqrt {d+e x}}{\sqrt {c d-b e}}\right )}{4 b^5 \sqrt {c} \sqrt {c d-b e}}-\frac {\sqrt {d+e x} \left (x \left (-b c (A e+B d)+2 A c^2 d+b^2 B e\right )+A b c d\right )}{2 b^2 c \left (b x+c x^2\right )^2}-\frac {\sqrt {d+e x} \left (b (c d-b e) \left (7 A b c e-12 A c^2 d-2 b^2 B e+6 b B c d\right )-3 c x (c d-b e) \left (-4 b c (A e+B d)+8 A c^2 d+b^2 B e\right )\right )}{4 b^4 c \left (b x+c x^2\right ) (c d-b e)} \]

[In]

Int[((A + B*x)*(d + e*x)^(3/2))/(b*x + c*x^2)^3,x]

[Out]

-1/2*(Sqrt[d + e*x]*(A*b*c*d + (2*A*c^2*d + b^2*B*e - b*c*(B*d + A*e))*x))/(b^2*c*(b*x + c*x^2)^2) - (Sqrt[d +
 e*x]*(b*(c*d - b*e)*(6*b*B*c*d - 12*A*c^2*d - 2*b^2*B*e + 7*A*b*c*e) - 3*c*(c*d - b*e)*(8*A*c^2*d + b^2*B*e -
 4*b*c*(B*d + A*e))*x))/(4*b^4*c*(c*d - b*e)*(b*x + c*x^2)) - (3*(16*A*c^2*d^2 + b^2*e*(4*B*d + A*e) - 4*b*c*d
*(2*B*d + 3*A*e))*ArcTanh[Sqrt[d + e*x]/Sqrt[d]])/(4*b^5*Sqrt[d]) + (3*(16*A*c^3*d^2 - b^3*B*e^2 - 4*b*c^2*d*(
2*B*d + 5*A*e) + b^2*c*e*(8*B*d + 5*A*e))*ArcTanh[(Sqrt[c]*Sqrt[d + e*x])/Sqrt[c*d - b*e]])/(4*b^5*Sqrt[c]*Sqr
t[c*d - b*e])

Rule 214

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-a/b, 2]/a)*ArcTanh[x/Rt[-a/b, 2]], x] /; FreeQ[{a, b},
x] && NegQ[a/b]

Rule 832

Int[((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_))*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_.), x_Symbol] :> Sim
p[(-(d + e*x)^(m - 1))*(a + b*x + c*x^2)^(p + 1)*((2*a*c*(e*f + d*g) - b*(c*d*f + a*e*g) - (2*c^2*d*f + b^2*e*
g - c*(b*e*f + b*d*g + 2*a*e*g))*x)/(c*(p + 1)*(b^2 - 4*a*c))), x] - Dist[1/(c*(p + 1)*(b^2 - 4*a*c)), Int[(d
+ e*x)^(m - 2)*(a + b*x + c*x^2)^(p + 1)*Simp[2*c^2*d^2*f*(2*p + 3) + b*e*g*(a*e*(m - 1) + b*d*(p + 2)) - c*(2
*a*e*(e*f*(m - 1) + d*g*m) + b*d*(d*g*(2*p + 3) - e*f*(m - 2*p - 4))) + e*(b^2*e*g*(m + p + 1) + 2*c^2*d*f*(m
+ 2*p + 2) - c*(2*a*e*g*m + b*(e*f + d*g)*(m + 2*p + 2)))*x, x], x], x] /; FreeQ[{a, b, c, d, e, f, g}, x] &&
NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0] && LtQ[p, -1] && GtQ[m, 1] && ((EqQ[m, 2] && EqQ[p, -3] &
& RationalQ[a, b, c, d, e, f, g]) ||  !ILtQ[m + 2*p + 3, 0])

Rule 836

Int[((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_))*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Simp
[(d + e*x)^(m + 1)*(f*(b*c*d - b^2*e + 2*a*c*e) - a*g*(2*c*d - b*e) + c*(f*(2*c*d - b*e) - g*(b*d - 2*a*e))*x)
*((a + b*x + c*x^2)^(p + 1)/((p + 1)*(b^2 - 4*a*c)*(c*d^2 - b*d*e + a*e^2))), x] + Dist[1/((p + 1)*(b^2 - 4*a*
c)*(c*d^2 - b*d*e + a*e^2)), Int[(d + e*x)^m*(a + b*x + c*x^2)^(p + 1)*Simp[f*(b*c*d*e*(2*p - m + 2) + b^2*e^2
*(p + m + 2) - 2*c^2*d^2*(2*p + 3) - 2*a*c*e^2*(m + 2*p + 3)) - g*(a*e*(b*e - 2*c*d*m + b*e*m) - b*d*(3*c*d -
b*e + 2*c*d*p - b*e*p)) + c*e*(g*(b*d - 2*a*e) - f*(2*c*d - b*e))*(m + 2*p + 4)*x, x], x], x] /; FreeQ[{a, b,
c, d, e, f, g, m}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0] && LtQ[p, -1] && (IntegerQ[m] ||
 IntegerQ[p] || IntegersQ[2*m, 2*p])

Rule 840

Int[((f_.) + (g_.)*(x_))/(Sqrt[(d_.) + (e_.)*(x_)]*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)), x_Symbol] :> Dist[2,
Subst[Int[(e*f - d*g + g*x^2)/(c*d^2 - b*d*e + a*e^2 - (2*c*d - b*e)*x^2 + c*x^4), x], x, Sqrt[d + e*x]], x] /
; FreeQ[{a, b, c, d, e, f, g}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0]

Rule 1180

Int[((d_) + (e_.)*(x_)^2)/((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[b^2 - 4*a*c, 2]}, Di
st[e/2 + (2*c*d - b*e)/(2*q), Int[1/(b/2 - q/2 + c*x^2), x], x] + Dist[e/2 - (2*c*d - b*e)/(2*q), Int[1/(b/2 +
 q/2 + c*x^2), x], x]] /; FreeQ[{a, b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - a*e^2, 0] && PosQ[b^
2 - 4*a*c]

Rubi steps \begin{align*} \text {integral}& = -\frac {\sqrt {d+e x} \left (A b c d+\left (2 A c^2 d+b^2 B e-b c (B d+A e)\right ) x\right )}{2 b^2 c \left (b x+c x^2\right )^2}+\frac {\int \frac {-\frac {1}{2} d \left (12 A c^2 d+2 b^2 B e-b c (6 B d+7 A e)\right )-\frac {1}{2} e \left (10 A c^2 d+b^2 B e-5 b c (B d+A e)\right ) x}{\sqrt {d+e x} \left (b x+c x^2\right )^2} \, dx}{2 b^2 c} \\ & = -\frac {\sqrt {d+e x} \left (A b c d+\left (2 A c^2 d+b^2 B e-b c (B d+A e)\right ) x\right )}{2 b^2 c \left (b x+c x^2\right )^2}-\frac {\sqrt {d+e x} \left (b (c d-b e) \left (6 b B c d-12 A c^2 d-2 b^2 B e+7 A b c e\right )-3 c (c d-b e) \left (8 A c^2 d+b^2 B e-4 b c (B d+A e)\right ) x\right )}{4 b^4 c (c d-b e) \left (b x+c x^2\right )}-\frac {\int \frac {-\frac {3}{4} c d (c d-b e) \left (16 A c^2 d^2+b^2 e (4 B d+A e)-4 b c d (2 B d+3 A e)\right )-\frac {3}{4} c d e (c d-b e) \left (8 A c^2 d+b^2 B e-4 b c (B d+A e)\right ) x}{\sqrt {d+e x} \left (b x+c x^2\right )} \, dx}{2 b^4 c d (c d-b e)} \\ & = -\frac {\sqrt {d+e x} \left (A b c d+\left (2 A c^2 d+b^2 B e-b c (B d+A e)\right ) x\right )}{2 b^2 c \left (b x+c x^2\right )^2}-\frac {\sqrt {d+e x} \left (b (c d-b e) \left (6 b B c d-12 A c^2 d-2 b^2 B e+7 A b c e\right )-3 c (c d-b e) \left (8 A c^2 d+b^2 B e-4 b c (B d+A e)\right ) x\right )}{4 b^4 c (c d-b e) \left (b x+c x^2\right )}-\frac {\text {Subst}\left (\int \frac {\frac {3}{4} c d^2 e (c d-b e) \left (8 A c^2 d+b^2 B e-4 b c (B d+A e)\right )-\frac {3}{4} c d e (c d-b e) \left (16 A c^2 d^2+b^2 e (4 B d+A e)-4 b c d (2 B d+3 A e)\right )-\frac {3}{4} c d e (c d-b e) \left (8 A c^2 d+b^2 B e-4 b c (B d+A e)\right ) x^2}{c d^2-b d e+(-2 c d+b e) x^2+c x^4} \, dx,x,\sqrt {d+e x}\right )}{b^4 c d (c d-b e)} \\ & = -\frac {\sqrt {d+e x} \left (A b c d+\left (2 A c^2 d+b^2 B e-b c (B d+A e)\right ) x\right )}{2 b^2 c \left (b x+c x^2\right )^2}-\frac {\sqrt {d+e x} \left (b (c d-b e) \left (6 b B c d-12 A c^2 d-2 b^2 B e+7 A b c e\right )-3 c (c d-b e) \left (8 A c^2 d+b^2 B e-4 b c (B d+A e)\right ) x\right )}{4 b^4 c (c d-b e) \left (b x+c x^2\right )}+\frac {\left (3 c \left (16 A c^2 d^2+b^2 e (4 B d+A e)-4 b c d (2 B d+3 A e)\right )\right ) \text {Subst}\left (\int \frac {1}{-\frac {b e}{2}+\frac {1}{2} (-2 c d+b e)+c x^2} \, dx,x,\sqrt {d+e x}\right )}{4 b^5}-\frac {\left (3 \left (16 A c^3 d^2-b^3 B e^2-4 b c^2 d (2 B d+5 A e)+b^2 c e (8 B d+5 A e)\right )\right ) \text {Subst}\left (\int \frac {1}{\frac {b e}{2}+\frac {1}{2} (-2 c d+b e)+c x^2} \, dx,x,\sqrt {d+e x}\right )}{4 b^5} \\ & = -\frac {\sqrt {d+e x} \left (A b c d+\left (2 A c^2 d+b^2 B e-b c (B d+A e)\right ) x\right )}{2 b^2 c \left (b x+c x^2\right )^2}-\frac {\sqrt {d+e x} \left (b (c d-b e) \left (6 b B c d-12 A c^2 d-2 b^2 B e+7 A b c e\right )-3 c (c d-b e) \left (8 A c^2 d+b^2 B e-4 b c (B d+A e)\right ) x\right )}{4 b^4 c (c d-b e) \left (b x+c x^2\right )}-\frac {3 \left (16 A c^2 d^2+b^2 e (4 B d+A e)-4 b c d (2 B d+3 A e)\right ) \tanh ^{-1}\left (\frac {\sqrt {d+e x}}{\sqrt {d}}\right )}{4 b^5 \sqrt {d}}+\frac {3 \left (16 A c^3 d^2-b^3 B e^2-4 b c^2 d (2 B d+5 A e)+b^2 c e (8 B d+5 A e)\right ) \tanh ^{-1}\left (\frac {\sqrt {c} \sqrt {d+e x}}{\sqrt {c d-b e}}\right )}{4 b^5 \sqrt {c} \sqrt {c d-b e}} \\ \end{align*}

Mathematica [A] (verified)

Time = 1.68 (sec) , antiderivative size = 279, normalized size of antiderivative = 0.81 \[ \int \frac {(A+B x) (d+e x)^{3/2}}{\left (b x+c x^2\right )^3} \, dx=-\frac {\frac {b \sqrt {d+e x} \left (b B x \left (12 c^2 d x^2+b^2 (4 d-5 e x)-3 b c x (-6 d+e x)\right )+A \left (-24 c^3 d x^3+12 b c^2 x^2 (-3 d+e x)+b^3 (2 d+5 e x)+b^2 c x (-8 d+19 e x)\right )\right )}{x^2 (b+c x)^2}-\frac {3 \left (-16 A c^3 d^2+b^3 B e^2+4 b c^2 d (2 B d+5 A e)-b^2 c e (8 B d+5 A e)\right ) \arctan \left (\frac {\sqrt {c} \sqrt {d+e x}}{\sqrt {-c d+b e}}\right )}{\sqrt {c} \sqrt {-c d+b e}}+\frac {3 \left (16 A c^2 d^2+b^2 e (4 B d+A e)-4 b c d (2 B d+3 A e)\right ) \text {arctanh}\left (\frac {\sqrt {d+e x}}{\sqrt {d}}\right )}{\sqrt {d}}}{4 b^5} \]

[In]

Integrate[((A + B*x)*(d + e*x)^(3/2))/(b*x + c*x^2)^3,x]

[Out]

-1/4*((b*Sqrt[d + e*x]*(b*B*x*(12*c^2*d*x^2 + b^2*(4*d - 5*e*x) - 3*b*c*x*(-6*d + e*x)) + A*(-24*c^3*d*x^3 + 1
2*b*c^2*x^2*(-3*d + e*x) + b^3*(2*d + 5*e*x) + b^2*c*x*(-8*d + 19*e*x))))/(x^2*(b + c*x)^2) - (3*(-16*A*c^3*d^
2 + b^3*B*e^2 + 4*b*c^2*d*(2*B*d + 5*A*e) - b^2*c*e*(8*B*d + 5*A*e))*ArcTan[(Sqrt[c]*Sqrt[d + e*x])/Sqrt[-(c*d
) + b*e]])/(Sqrt[c]*Sqrt[-(c*d) + b*e]) + (3*(16*A*c^2*d^2 + b^2*e*(4*B*d + A*e) - 4*b*c*d*(2*B*d + 3*A*e))*Ar
cTanh[Sqrt[d + e*x]/Sqrt[d]])/Sqrt[d])/b^5

Maple [A] (verified)

Time = 0.57 (sec) , antiderivative size = 314, normalized size of antiderivative = 0.91

method result size
pseudoelliptic \(-\frac {12 \left (x^{2} \left (c x +b \right )^{2} \left (\frac {5 \left (A c -\frac {B b}{5}\right ) e^{2} b^{2} \sqrt {d}}{16}+c \left (\frac {b^{2} B e}{2}-\frac {5 b \left (A e +\frac {2 B d}{5}\right ) c}{4}+A \,c^{2} d \right ) d^{\frac {3}{2}}\right ) \arctan \left (\frac {c \sqrt {e x +d}}{\sqrt {\left (b e -c d \right ) c}}\right )+\frac {\sqrt {\left (b e -c d \right ) c}\, \left (\frac {3 x^{2} \left (c x +b \right )^{2} \left (\left (A \,e^{2}+4 B d e \right ) b^{2}+4 \left (-3 A d e -2 B \,d^{2}\right ) b c +16 A \,c^{2} d^{2}\right ) \operatorname {arctanh}\left (\frac {\sqrt {e x +d}}{\sqrt {d}}\right )}{2}+\sqrt {e x +d}\, \left (\frac {5 \left (\frac {3 \left (4 A \,c^{2}-B b c \right ) x^{2}}{5}+\left (\frac {19}{5} A b c -B \,b^{2}\right ) x +A \,b^{2}\right ) x e b \sqrt {d}}{2}+\left (6 \left (-2 A \,c^{3}+B b \,c^{2}\right ) x^{3}+9 \left (-2 A b \,c^{2}+B \,b^{2} c \right ) x^{2}+2 \left (-2 A \,b^{2} c +B \,b^{3}\right ) x +A \,b^{3}\right ) d^{\frac {3}{2}}\right ) b \right )}{24}\right )}{\sqrt {\left (b e -c d \right ) c}\, \sqrt {d}\, \left (c x +b \right )^{2} b^{5} x^{2}}\) \(314\)
risch \(-\frac {\sqrt {e x +d}\, \left (5 A b e x -12 A c d x +4 B b d x +2 A b d \right )}{4 b^{4} x^{2}}-\frac {e \left (\frac {\frac {8 \left (\left (\frac {7}{8} A \,b^{2} c^{2} e^{2}-\frac {3}{2} A b \,c^{3} d e -\frac {3}{8} B \,b^{3} c \,e^{2}+B \,b^{2} c^{2} d e \right ) \left (e x +d \right )^{\frac {3}{2}}+\frac {e b \left (9 A \,b^{2} c \,e^{2}-21 A b \,c^{2} d e +12 A \,c^{3} d^{2}-5 b^{3} B \,e^{2}+13 B \,b^{2} c d e -8 B b \,c^{2} d^{2}\right ) \sqrt {e x +d}}{8}\right )}{\left (c \left (e x +d \right )+b e -c d \right )^{2}}+\frac {3 \left (5 A \,b^{2} c \,e^{2}-20 A b \,c^{2} d e +16 A \,c^{3} d^{2}-b^{3} B \,e^{2}+8 B \,b^{2} c d e -8 B b \,c^{2} d^{2}\right ) \arctan \left (\frac {c \sqrt {e x +d}}{\sqrt {\left (b e -c d \right ) c}}\right )}{\sqrt {\left (b e -c d \right ) c}}}{b e}-\frac {\left (-3 A \,b^{2} e^{2}+36 A b c d e -48 A \,c^{2} d^{2}-12 B \,b^{2} d e +24 c \,d^{2} B b \right ) \operatorname {arctanh}\left (\frac {\sqrt {e x +d}}{\sqrt {d}}\right )}{b e \sqrt {d}}\right )}{4 b^{4}}\) \(352\)
derivativedivides \(2 e^{4} \left (-\frac {\frac {\left (\frac {7}{8} A \,b^{2} c^{2} e^{2}-\frac {3}{2} A b \,c^{3} d e -\frac {3}{8} B \,b^{3} c \,e^{2}+B \,b^{2} c^{2} d e \right ) \left (e x +d \right )^{\frac {3}{2}}+\frac {e b \left (9 A \,b^{2} c \,e^{2}-21 A b \,c^{2} d e +12 A \,c^{3} d^{2}-5 b^{3} B \,e^{2}+13 B \,b^{2} c d e -8 B b \,c^{2} d^{2}\right ) \sqrt {e x +d}}{8}}{\left (c \left (e x +d \right )+b e -c d \right )^{2}}+\frac {3 \left (5 A \,b^{2} c \,e^{2}-20 A b \,c^{2} d e +16 A \,c^{3} d^{2}-b^{3} B \,e^{2}+8 B \,b^{2} c d e -8 B b \,c^{2} d^{2}\right ) \arctan \left (\frac {c \sqrt {e x +d}}{\sqrt {\left (b e -c d \right ) c}}\right )}{8 \sqrt {\left (b e -c d \right ) c}}}{e^{4} b^{5}}-\frac {\frac {\frac {e b \left (5 A b e -12 A c d +4 B b d \right ) \left (e x +d \right )^{\frac {3}{2}}}{8}+\left (\frac {3}{2} A b c \,d^{2} e -\frac {1}{2} B \,b^{2} d^{2} e -\frac {3}{8} A \,b^{2} d \,e^{2}\right ) \sqrt {e x +d}}{e^{2} x^{2}}+\frac {3 \left (A \,b^{2} e^{2}-12 A b c d e +16 A \,c^{2} d^{2}+4 B \,b^{2} d e -8 c \,d^{2} B b \right ) \operatorname {arctanh}\left (\frac {\sqrt {e x +d}}{\sqrt {d}}\right )}{8 \sqrt {d}}}{e^{4} b^{5}}\right )\) \(386\)
default \(2 e^{4} \left (-\frac {\frac {\left (\frac {7}{8} A \,b^{2} c^{2} e^{2}-\frac {3}{2} A b \,c^{3} d e -\frac {3}{8} B \,b^{3} c \,e^{2}+B \,b^{2} c^{2} d e \right ) \left (e x +d \right )^{\frac {3}{2}}+\frac {e b \left (9 A \,b^{2} c \,e^{2}-21 A b \,c^{2} d e +12 A \,c^{3} d^{2}-5 b^{3} B \,e^{2}+13 B \,b^{2} c d e -8 B b \,c^{2} d^{2}\right ) \sqrt {e x +d}}{8}}{\left (c \left (e x +d \right )+b e -c d \right )^{2}}+\frac {3 \left (5 A \,b^{2} c \,e^{2}-20 A b \,c^{2} d e +16 A \,c^{3} d^{2}-b^{3} B \,e^{2}+8 B \,b^{2} c d e -8 B b \,c^{2} d^{2}\right ) \arctan \left (\frac {c \sqrt {e x +d}}{\sqrt {\left (b e -c d \right ) c}}\right )}{8 \sqrt {\left (b e -c d \right ) c}}}{e^{4} b^{5}}-\frac {\frac {\frac {e b \left (5 A b e -12 A c d +4 B b d \right ) \left (e x +d \right )^{\frac {3}{2}}}{8}+\left (\frac {3}{2} A b c \,d^{2} e -\frac {1}{2} B \,b^{2} d^{2} e -\frac {3}{8} A \,b^{2} d \,e^{2}\right ) \sqrt {e x +d}}{e^{2} x^{2}}+\frac {3 \left (A \,b^{2} e^{2}-12 A b c d e +16 A \,c^{2} d^{2}+4 B \,b^{2} d e -8 c \,d^{2} B b \right ) \operatorname {arctanh}\left (\frac {\sqrt {e x +d}}{\sqrt {d}}\right )}{8 \sqrt {d}}}{e^{4} b^{5}}\right )\) \(386\)

[In]

int((B*x+A)*(e*x+d)^(3/2)/(c*x^2+b*x)^3,x,method=_RETURNVERBOSE)

[Out]

-12*(x^2*(c*x+b)^2*(5/16*(A*c-1/5*B*b)*e^2*b^2*d^(1/2)+c*(1/2*b^2*B*e-5/4*b*(A*e+2/5*B*d)*c+A*c^2*d)*d^(3/2))*
arctan(c*(e*x+d)^(1/2)/((b*e-c*d)*c)^(1/2))+1/24*((b*e-c*d)*c)^(1/2)*(3/2*x^2*(c*x+b)^2*((A*e^2+4*B*d*e)*b^2+4
*(-3*A*d*e-2*B*d^2)*b*c+16*A*c^2*d^2)*arctanh((e*x+d)^(1/2)/d^(1/2))+(e*x+d)^(1/2)*(5/2*(3/5*(4*A*c^2-B*b*c)*x
^2+(19/5*A*b*c-B*b^2)*x+A*b^2)*x*e*b*d^(1/2)+(6*(-2*A*c^3+B*b*c^2)*x^3+9*(-2*A*b*c^2+B*b^2*c)*x^2+2*(-2*A*b^2*
c+B*b^3)*x+A*b^3)*d^(3/2))*b))/((b*e-c*d)*c)^(1/2)/d^(1/2)/(c*x+b)^2/b^5/x^2

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 859 vs. \(2 (319) = 638\).

Time = 1.43 (sec) , antiderivative size = 3471, normalized size of antiderivative = 10.03 \[ \int \frac {(A+B x) (d+e x)^{3/2}}{\left (b x+c x^2\right )^3} \, dx=\text {Too large to display} \]

[In]

integrate((B*x+A)*(e*x+d)^(3/2)/(c*x^2+b*x)^3,x, algorithm="fricas")

[Out]

[1/8*(3*((8*(B*b*c^4 - 2*A*c^5)*d^3 - 4*(2*B*b^2*c^3 - 5*A*b*c^4)*d^2*e + (B*b^3*c^2 - 5*A*b^2*c^3)*d*e^2)*x^4
 + 2*(8*(B*b^2*c^3 - 2*A*b*c^4)*d^3 - 4*(2*B*b^3*c^2 - 5*A*b^2*c^3)*d^2*e + (B*b^4*c - 5*A*b^3*c^2)*d*e^2)*x^3
 + (8*(B*b^3*c^2 - 2*A*b^2*c^3)*d^3 - 4*(2*B*b^4*c - 5*A*b^3*c^2)*d^2*e + (B*b^5 - 5*A*b^4*c)*d*e^2)*x^2)*sqrt
(c^2*d - b*c*e)*log((c*e*x + 2*c*d - b*e - 2*sqrt(c^2*d - b*c*e)*sqrt(e*x + d))/(c*x + b)) - 3*((A*b^3*c^3*e^3
 + 8*(B*b*c^5 - 2*A*c^6)*d^3 - 4*(3*B*b^2*c^4 - 7*A*b*c^5)*d^2*e + (4*B*b^3*c^3 - 13*A*b^2*c^4)*d*e^2)*x^4 + 2
*(A*b^4*c^2*e^3 + 8*(B*b^2*c^4 - 2*A*b*c^5)*d^3 - 4*(3*B*b^3*c^3 - 7*A*b^2*c^4)*d^2*e + (4*B*b^4*c^2 - 13*A*b^
3*c^3)*d*e^2)*x^3 + (A*b^5*c*e^3 + 8*(B*b^3*c^3 - 2*A*b^2*c^4)*d^3 - 4*(3*B*b^4*c^2 - 7*A*b^3*c^3)*d^2*e + (4*
B*b^5*c - 13*A*b^4*c^2)*d*e^2)*x^2)*sqrt(d)*log((e*x - 2*sqrt(e*x + d)*sqrt(d) + 2*d)/x) - 2*(2*A*b^4*c^2*d^3
- 2*A*b^5*c*d^2*e + 3*(4*(B*b^2*c^4 - 2*A*b*c^5)*d^3 - (5*B*b^3*c^3 - 12*A*b^2*c^4)*d^2*e + (B*b^4*c^2 - 4*A*b
^3*c^3)*d*e^2)*x^3 + (18*(B*b^3*c^3 - 2*A*b^2*c^4)*d^3 - (23*B*b^4*c^2 - 55*A*b^3*c^3)*d^2*e + (5*B*b^5*c - 19
*A*b^4*c^2)*d*e^2)*x^2 - (5*A*b^5*c*d*e^2 - 4*(B*b^4*c^2 - 2*A*b^3*c^3)*d^3 + (4*B*b^5*c - 13*A*b^4*c^2)*d^2*e
)*x)*sqrt(e*x + d))/((b^5*c^4*d^2 - b^6*c^3*d*e)*x^4 + 2*(b^6*c^3*d^2 - b^7*c^2*d*e)*x^3 + (b^7*c^2*d^2 - b^8*
c*d*e)*x^2), 1/8*(6*((8*(B*b*c^4 - 2*A*c^5)*d^3 - 4*(2*B*b^2*c^3 - 5*A*b*c^4)*d^2*e + (B*b^3*c^2 - 5*A*b^2*c^3
)*d*e^2)*x^4 + 2*(8*(B*b^2*c^3 - 2*A*b*c^4)*d^3 - 4*(2*B*b^3*c^2 - 5*A*b^2*c^3)*d^2*e + (B*b^4*c - 5*A*b^3*c^2
)*d*e^2)*x^3 + (8*(B*b^3*c^2 - 2*A*b^2*c^3)*d^3 - 4*(2*B*b^4*c - 5*A*b^3*c^2)*d^2*e + (B*b^5 - 5*A*b^4*c)*d*e^
2)*x^2)*sqrt(-c^2*d + b*c*e)*arctan(sqrt(-c^2*d + b*c*e)*sqrt(e*x + d)/(c*e*x + c*d)) - 3*((A*b^3*c^3*e^3 + 8*
(B*b*c^5 - 2*A*c^6)*d^3 - 4*(3*B*b^2*c^4 - 7*A*b*c^5)*d^2*e + (4*B*b^3*c^3 - 13*A*b^2*c^4)*d*e^2)*x^4 + 2*(A*b
^4*c^2*e^3 + 8*(B*b^2*c^4 - 2*A*b*c^5)*d^3 - 4*(3*B*b^3*c^3 - 7*A*b^2*c^4)*d^2*e + (4*B*b^4*c^2 - 13*A*b^3*c^3
)*d*e^2)*x^3 + (A*b^5*c*e^3 + 8*(B*b^3*c^3 - 2*A*b^2*c^4)*d^3 - 4*(3*B*b^4*c^2 - 7*A*b^3*c^3)*d^2*e + (4*B*b^5
*c - 13*A*b^4*c^2)*d*e^2)*x^2)*sqrt(d)*log((e*x - 2*sqrt(e*x + d)*sqrt(d) + 2*d)/x) - 2*(2*A*b^4*c^2*d^3 - 2*A
*b^5*c*d^2*e + 3*(4*(B*b^2*c^4 - 2*A*b*c^5)*d^3 - (5*B*b^3*c^3 - 12*A*b^2*c^4)*d^2*e + (B*b^4*c^2 - 4*A*b^3*c^
3)*d*e^2)*x^3 + (18*(B*b^3*c^3 - 2*A*b^2*c^4)*d^3 - (23*B*b^4*c^2 - 55*A*b^3*c^3)*d^2*e + (5*B*b^5*c - 19*A*b^
4*c^2)*d*e^2)*x^2 - (5*A*b^5*c*d*e^2 - 4*(B*b^4*c^2 - 2*A*b^3*c^3)*d^3 + (4*B*b^5*c - 13*A*b^4*c^2)*d^2*e)*x)*
sqrt(e*x + d))/((b^5*c^4*d^2 - b^6*c^3*d*e)*x^4 + 2*(b^6*c^3*d^2 - b^7*c^2*d*e)*x^3 + (b^7*c^2*d^2 - b^8*c*d*e
)*x^2), -1/8*(6*((A*b^3*c^3*e^3 + 8*(B*b*c^5 - 2*A*c^6)*d^3 - 4*(3*B*b^2*c^4 - 7*A*b*c^5)*d^2*e + (4*B*b^3*c^3
 - 13*A*b^2*c^4)*d*e^2)*x^4 + 2*(A*b^4*c^2*e^3 + 8*(B*b^2*c^4 - 2*A*b*c^5)*d^3 - 4*(3*B*b^3*c^3 - 7*A*b^2*c^4)
*d^2*e + (4*B*b^4*c^2 - 13*A*b^3*c^3)*d*e^2)*x^3 + (A*b^5*c*e^3 + 8*(B*b^3*c^3 - 2*A*b^2*c^4)*d^3 - 4*(3*B*b^4
*c^2 - 7*A*b^3*c^3)*d^2*e + (4*B*b^5*c - 13*A*b^4*c^2)*d*e^2)*x^2)*sqrt(-d)*arctan(sqrt(e*x + d)*sqrt(-d)/d) -
 3*((8*(B*b*c^4 - 2*A*c^5)*d^3 - 4*(2*B*b^2*c^3 - 5*A*b*c^4)*d^2*e + (B*b^3*c^2 - 5*A*b^2*c^3)*d*e^2)*x^4 + 2*
(8*(B*b^2*c^3 - 2*A*b*c^4)*d^3 - 4*(2*B*b^3*c^2 - 5*A*b^2*c^3)*d^2*e + (B*b^4*c - 5*A*b^3*c^2)*d*e^2)*x^3 + (8
*(B*b^3*c^2 - 2*A*b^2*c^3)*d^3 - 4*(2*B*b^4*c - 5*A*b^3*c^2)*d^2*e + (B*b^5 - 5*A*b^4*c)*d*e^2)*x^2)*sqrt(c^2*
d - b*c*e)*log((c*e*x + 2*c*d - b*e - 2*sqrt(c^2*d - b*c*e)*sqrt(e*x + d))/(c*x + b)) + 2*(2*A*b^4*c^2*d^3 - 2
*A*b^5*c*d^2*e + 3*(4*(B*b^2*c^4 - 2*A*b*c^5)*d^3 - (5*B*b^3*c^3 - 12*A*b^2*c^4)*d^2*e + (B*b^4*c^2 - 4*A*b^3*
c^3)*d*e^2)*x^3 + (18*(B*b^3*c^3 - 2*A*b^2*c^4)*d^3 - (23*B*b^4*c^2 - 55*A*b^3*c^3)*d^2*e + (5*B*b^5*c - 19*A*
b^4*c^2)*d*e^2)*x^2 - (5*A*b^5*c*d*e^2 - 4*(B*b^4*c^2 - 2*A*b^3*c^3)*d^3 + (4*B*b^5*c - 13*A*b^4*c^2)*d^2*e)*x
)*sqrt(e*x + d))/((b^5*c^4*d^2 - b^6*c^3*d*e)*x^4 + 2*(b^6*c^3*d^2 - b^7*c^2*d*e)*x^3 + (b^7*c^2*d^2 - b^8*c*d
*e)*x^2), 1/4*(3*((8*(B*b*c^4 - 2*A*c^5)*d^3 - 4*(2*B*b^2*c^3 - 5*A*b*c^4)*d^2*e + (B*b^3*c^2 - 5*A*b^2*c^3)*d
*e^2)*x^4 + 2*(8*(B*b^2*c^3 - 2*A*b*c^4)*d^3 - 4*(2*B*b^3*c^2 - 5*A*b^2*c^3)*d^2*e + (B*b^4*c - 5*A*b^3*c^2)*d
*e^2)*x^3 + (8*(B*b^3*c^2 - 2*A*b^2*c^3)*d^3 - 4*(2*B*b^4*c - 5*A*b^3*c^2)*d^2*e + (B*b^5 - 5*A*b^4*c)*d*e^2)*
x^2)*sqrt(-c^2*d + b*c*e)*arctan(sqrt(-c^2*d + b*c*e)*sqrt(e*x + d)/(c*e*x + c*d)) - 3*((A*b^3*c^3*e^3 + 8*(B*
b*c^5 - 2*A*c^6)*d^3 - 4*(3*B*b^2*c^4 - 7*A*b*c^5)*d^2*e + (4*B*b^3*c^3 - 13*A*b^2*c^4)*d*e^2)*x^4 + 2*(A*b^4*
c^2*e^3 + 8*(B*b^2*c^4 - 2*A*b*c^5)*d^3 - 4*(3*B*b^3*c^3 - 7*A*b^2*c^4)*d^2*e + (4*B*b^4*c^2 - 13*A*b^3*c^3)*d
*e^2)*x^3 + (A*b^5*c*e^3 + 8*(B*b^3*c^3 - 2*A*b^2*c^4)*d^3 - 4*(3*B*b^4*c^2 - 7*A*b^3*c^3)*d^2*e + (4*B*b^5*c
- 13*A*b^4*c^2)*d*e^2)*x^2)*sqrt(-d)*arctan(sqrt(e*x + d)*sqrt(-d)/d) - (2*A*b^4*c^2*d^3 - 2*A*b^5*c*d^2*e + 3
*(4*(B*b^2*c^4 - 2*A*b*c^5)*d^3 - (5*B*b^3*c^3 - 12*A*b^2*c^4)*d^2*e + (B*b^4*c^2 - 4*A*b^3*c^3)*d*e^2)*x^3 +
(18*(B*b^3*c^3 - 2*A*b^2*c^4)*d^3 - (23*B*b^4*c^2 - 55*A*b^3*c^3)*d^2*e + (5*B*b^5*c - 19*A*b^4*c^2)*d*e^2)*x^
2 - (5*A*b^5*c*d*e^2 - 4*(B*b^4*c^2 - 2*A*b^3*c^3)*d^3 + (4*B*b^5*c - 13*A*b^4*c^2)*d^2*e)*x)*sqrt(e*x + d))/(
(b^5*c^4*d^2 - b^6*c^3*d*e)*x^4 + 2*(b^6*c^3*d^2 - b^7*c^2*d*e)*x^3 + (b^7*c^2*d^2 - b^8*c*d*e)*x^2)]

Sympy [F(-1)]

Timed out. \[ \int \frac {(A+B x) (d+e x)^{3/2}}{\left (b x+c x^2\right )^3} \, dx=\text {Timed out} \]

[In]

integrate((B*x+A)*(e*x+d)**(3/2)/(c*x**2+b*x)**3,x)

[Out]

Timed out

Maxima [F(-2)]

Exception generated. \[ \int \frac {(A+B x) (d+e x)^{3/2}}{\left (b x+c x^2\right )^3} \, dx=\text {Exception raised: ValueError} \]

[In]

integrate((B*x+A)*(e*x+d)^(3/2)/(c*x^2+b*x)^3,x, algorithm="maxima")

[Out]

Exception raised: ValueError >> Computation failed since Maxima requested additional constraints; using the 'a
ssume' command before evaluation *may* help (example of legal syntax is 'assume(b*e-c*d>0)', see `assume?` for
 more detail

Giac [A] (verification not implemented)

none

Time = 0.29 (sec) , antiderivative size = 635, normalized size of antiderivative = 1.84 \[ \int \frac {(A+B x) (d+e x)^{3/2}}{\left (b x+c x^2\right )^3} \, dx=\frac {3 \, {\left (8 \, B b c^{2} d^{2} - 16 \, A c^{3} d^{2} - 8 \, B b^{2} c d e + 20 \, A b c^{2} d e + B b^{3} e^{2} - 5 \, A b^{2} c e^{2}\right )} \arctan \left (\frac {\sqrt {e x + d} c}{\sqrt {-c^{2} d + b c e}}\right )}{4 \, \sqrt {-c^{2} d + b c e} b^{5}} - \frac {3 \, {\left (8 \, B b c d^{2} - 16 \, A c^{2} d^{2} - 4 \, B b^{2} d e + 12 \, A b c d e - A b^{2} e^{2}\right )} \arctan \left (\frac {\sqrt {e x + d}}{\sqrt {-d}}\right )}{4 \, b^{5} \sqrt {-d}} - \frac {12 \, {\left (e x + d\right )}^{\frac {7}{2}} B b c^{2} d e - 24 \, {\left (e x + d\right )}^{\frac {7}{2}} A c^{3} d e - 36 \, {\left (e x + d\right )}^{\frac {5}{2}} B b c^{2} d^{2} e + 72 \, {\left (e x + d\right )}^{\frac {5}{2}} A c^{3} d^{2} e + 36 \, {\left (e x + d\right )}^{\frac {3}{2}} B b c^{2} d^{3} e - 72 \, {\left (e x + d\right )}^{\frac {3}{2}} A c^{3} d^{3} e - 12 \, \sqrt {e x + d} B b c^{2} d^{4} e + 24 \, \sqrt {e x + d} A c^{3} d^{4} e - 3 \, {\left (e x + d\right )}^{\frac {7}{2}} B b^{2} c e^{2} + 12 \, {\left (e x + d\right )}^{\frac {7}{2}} A b c^{2} e^{2} + 27 \, {\left (e x + d\right )}^{\frac {5}{2}} B b^{2} c d e^{2} - 72 \, {\left (e x + d\right )}^{\frac {5}{2}} A b c^{2} d e^{2} - 45 \, {\left (e x + d\right )}^{\frac {3}{2}} B b^{2} c d^{2} e^{2} + 108 \, {\left (e x + d\right )}^{\frac {3}{2}} A b c^{2} d^{2} e^{2} + 21 \, \sqrt {e x + d} B b^{2} c d^{3} e^{2} - 48 \, \sqrt {e x + d} A b c^{2} d^{3} e^{2} - 5 \, {\left (e x + d\right )}^{\frac {5}{2}} B b^{3} e^{3} + 19 \, {\left (e x + d\right )}^{\frac {5}{2}} A b^{2} c e^{3} + 14 \, {\left (e x + d\right )}^{\frac {3}{2}} B b^{3} d e^{3} - 46 \, {\left (e x + d\right )}^{\frac {3}{2}} A b^{2} c d e^{3} - 9 \, \sqrt {e x + d} B b^{3} d^{2} e^{3} + 27 \, \sqrt {e x + d} A b^{2} c d^{2} e^{3} + 5 \, {\left (e x + d\right )}^{\frac {3}{2}} A b^{3} e^{4} - 3 \, \sqrt {e x + d} A b^{3} d e^{4}}{4 \, {\left ({\left (e x + d\right )}^{2} c - 2 \, {\left (e x + d\right )} c d + c d^{2} + {\left (e x + d\right )} b e - b d e\right )}^{2} b^{4}} \]

[In]

integrate((B*x+A)*(e*x+d)^(3/2)/(c*x^2+b*x)^3,x, algorithm="giac")

[Out]

3/4*(8*B*b*c^2*d^2 - 16*A*c^3*d^2 - 8*B*b^2*c*d*e + 20*A*b*c^2*d*e + B*b^3*e^2 - 5*A*b^2*c*e^2)*arctan(sqrt(e*
x + d)*c/sqrt(-c^2*d + b*c*e))/(sqrt(-c^2*d + b*c*e)*b^5) - 3/4*(8*B*b*c*d^2 - 16*A*c^2*d^2 - 4*B*b^2*d*e + 12
*A*b*c*d*e - A*b^2*e^2)*arctan(sqrt(e*x + d)/sqrt(-d))/(b^5*sqrt(-d)) - 1/4*(12*(e*x + d)^(7/2)*B*b*c^2*d*e -
24*(e*x + d)^(7/2)*A*c^3*d*e - 36*(e*x + d)^(5/2)*B*b*c^2*d^2*e + 72*(e*x + d)^(5/2)*A*c^3*d^2*e + 36*(e*x + d
)^(3/2)*B*b*c^2*d^3*e - 72*(e*x + d)^(3/2)*A*c^3*d^3*e - 12*sqrt(e*x + d)*B*b*c^2*d^4*e + 24*sqrt(e*x + d)*A*c
^3*d^4*e - 3*(e*x + d)^(7/2)*B*b^2*c*e^2 + 12*(e*x + d)^(7/2)*A*b*c^2*e^2 + 27*(e*x + d)^(5/2)*B*b^2*c*d*e^2 -
 72*(e*x + d)^(5/2)*A*b*c^2*d*e^2 - 45*(e*x + d)^(3/2)*B*b^2*c*d^2*e^2 + 108*(e*x + d)^(3/2)*A*b*c^2*d^2*e^2 +
 21*sqrt(e*x + d)*B*b^2*c*d^3*e^2 - 48*sqrt(e*x + d)*A*b*c^2*d^3*e^2 - 5*(e*x + d)^(5/2)*B*b^3*e^3 + 19*(e*x +
 d)^(5/2)*A*b^2*c*e^3 + 14*(e*x + d)^(3/2)*B*b^3*d*e^3 - 46*(e*x + d)^(3/2)*A*b^2*c*d*e^3 - 9*sqrt(e*x + d)*B*
b^3*d^2*e^3 + 27*sqrt(e*x + d)*A*b^2*c*d^2*e^3 + 5*(e*x + d)^(3/2)*A*b^3*e^4 - 3*sqrt(e*x + d)*A*b^3*d*e^4)/((
(e*x + d)^2*c - 2*(e*x + d)*c*d + c*d^2 + (e*x + d)*b*e - b*d*e)^2*b^4)

Mupad [B] (verification not implemented)

Time = 13.09 (sec) , antiderivative size = 5796, normalized size of antiderivative = 16.75 \[ \int \frac {(A+B x) (d+e x)^{3/2}}{\left (b x+c x^2\right )^3} \, dx=\text {Too large to display} \]

[In]

int(((A + B*x)*(d + e*x)^(3/2))/(b*x + c*x^2)^3,x)

[Out]

((3*(d + e*x)^(1/2)*(A*b^3*d*e^4 - 8*A*c^3*d^4*e + 3*B*b^3*d^2*e^3 + 16*A*b*c^2*d^3*e^2 - 9*A*b^2*c*d^2*e^3 -
7*B*b^2*c*d^3*e^2 + 4*B*b*c^2*d^4*e))/(4*b^4) - ((d + e*x)^(3/2)*(5*A*b^3*e^4 - 72*A*c^3*d^3*e + 14*B*b^3*d*e^
3 + 108*A*b*c^2*d^2*e^2 - 45*B*b^2*c*d^2*e^2 - 46*A*b^2*c*d*e^3 + 36*B*b*c^2*d^3*e))/(4*b^4) + ((d + e*x)^(5/2
)*(5*B*b^3*e^3 - 19*A*b^2*c*e^3 - 72*A*c^3*d^2*e + 72*A*b*c^2*d*e^2 + 36*B*b*c^2*d^2*e - 27*B*b^2*c*d*e^2))/(4
*b^4) + (3*c*(d + e*x)^(7/2)*(B*b^2*e^2 - 4*A*b*c*e^2 + 8*A*c^2*d*e - 4*B*b*c*d*e))/(4*b^4))/(c^2*(d + e*x)^4
- (d + e*x)*(4*c^2*d^3 + 2*b^2*d*e^2 - 6*b*c*d^2*e) - (4*c^2*d - 2*b*c*e)*(d + e*x)^3 + (d + e*x)^2*(b^2*e^2 +
 6*c^2*d^2 - 6*b*c*d*e) + c^2*d^4 + b^2*d^2*e^2 - 2*b*c*d^3*e) + (atan((((-c*(b*e - c*d))^(1/2)*(((d + e*x)^(1
/2)*(234*A^2*b^4*c^3*e^6 + 4608*A^2*c^7*d^4*e^2 + 9*B^2*b^6*c*e^6 + 6624*A^2*b^2*c^5*d^2*e^4 + 1152*B^2*b^2*c^
5*d^4*e^2 - 1728*B^2*b^3*c^4*d^3*e^3 + 864*B^2*b^4*c^3*d^2*e^4 - 90*A*B*b^5*c^2*e^6 - 9216*A^2*b*c^6*d^3*e^3 -
 2016*A^2*b^3*c^4*d*e^5 - 144*B^2*b^5*c^2*d*e^5 - 4608*A*B*b*c^6*d^4*e^2 + 1152*A*B*b^4*c^3*d*e^5 + 8064*A*B*b
^2*c^5*d^3*e^3 - 4896*A*B*b^3*c^4*d^2*e^4))/(8*b^8) - (3*(-c*(b*e - c*d))^(1/2)*((3*A*b^12*c^2*e^5 - 24*A*b^11
*c^3*d*e^4 + 9*B*b^12*c^2*d*e^4 + 24*A*b^10*c^4*d^2*e^3 - 12*B*b^11*c^3*d^2*e^3)/b^12 - (3*(64*b^11*c^2*e^3 -
128*b^10*c^3*d*e^2)*(-c*(b*e - c*d))^(1/2)*(d + e*x)^(1/2)*(16*A*c^3*d^2 - B*b^3*e^2 + 5*A*b^2*c*e^2 - 8*B*b*c
^2*d^2 - 20*A*b*c^2*d*e + 8*B*b^2*c*d*e))/(64*b^8*(b^5*c^2*d - b^6*c*e)))*(16*A*c^3*d^2 - B*b^3*e^2 + 5*A*b^2*
c*e^2 - 8*B*b*c^2*d^2 - 20*A*b*c^2*d*e + 8*B*b^2*c*d*e))/(8*(b^5*c^2*d - b^6*c*e)))*(16*A*c^3*d^2 - B*b^3*e^2
+ 5*A*b^2*c*e^2 - 8*B*b*c^2*d^2 - 20*A*b*c^2*d*e + 8*B*b^2*c*d*e)*3i)/(8*(b^5*c^2*d - b^6*c*e)) + ((-c*(b*e -
c*d))^(1/2)*(((d + e*x)^(1/2)*(234*A^2*b^4*c^3*e^6 + 4608*A^2*c^7*d^4*e^2 + 9*B^2*b^6*c*e^6 + 6624*A^2*b^2*c^5
*d^2*e^4 + 1152*B^2*b^2*c^5*d^4*e^2 - 1728*B^2*b^3*c^4*d^3*e^3 + 864*B^2*b^4*c^3*d^2*e^4 - 90*A*B*b^5*c^2*e^6
- 9216*A^2*b*c^6*d^3*e^3 - 2016*A^2*b^3*c^4*d*e^5 - 144*B^2*b^5*c^2*d*e^5 - 4608*A*B*b*c^6*d^4*e^2 + 1152*A*B*
b^4*c^3*d*e^5 + 8064*A*B*b^2*c^5*d^3*e^3 - 4896*A*B*b^3*c^4*d^2*e^4))/(8*b^8) + (3*(-c*(b*e - c*d))^(1/2)*((3*
A*b^12*c^2*e^5 - 24*A*b^11*c^3*d*e^4 + 9*B*b^12*c^2*d*e^4 + 24*A*b^10*c^4*d^2*e^3 - 12*B*b^11*c^3*d^2*e^3)/b^1
2 + (3*(64*b^11*c^2*e^3 - 128*b^10*c^3*d*e^2)*(-c*(b*e - c*d))^(1/2)*(d + e*x)^(1/2)*(16*A*c^3*d^2 - B*b^3*e^2
 + 5*A*b^2*c*e^2 - 8*B*b*c^2*d^2 - 20*A*b*c^2*d*e + 8*B*b^2*c*d*e))/(64*b^8*(b^5*c^2*d - b^6*c*e)))*(16*A*c^3*
d^2 - B*b^3*e^2 + 5*A*b^2*c*e^2 - 8*B*b*c^2*d^2 - 20*A*b*c^2*d*e + 8*B*b^2*c*d*e))/(8*(b^5*c^2*d - b^6*c*e)))*
(16*A*c^3*d^2 - B*b^3*e^2 + 5*A*b^2*c*e^2 - 8*B*b*c^2*d^2 - 20*A*b*c^2*d*e + 8*B*b^2*c*d*e)*3i)/(8*(b^5*c^2*d
- b^6*c*e)))/(((135*A^3*b^5*c^3*e^8)/8 - 1728*A^3*c^8*d^5*e^3 - 3996*A^3*b^2*c^6*d^3*e^5 + 1674*A^3*b^3*c^5*d^
2*e^6 + 216*B^3*b^3*c^5*d^5*e^3 - 378*B^3*b^4*c^4*d^4*e^4 + 216*B^3*b^5*c^3*d^3*e^5 - (189*B^3*b^6*c^2*d^2*e^6
)/4 + (27*A*B^2*b^7*c*e^8)/32 + (27*B^3*b^7*c*d*e^7)/8 - (243*A^2*B*b^6*c^2*e^8)/32 + 4320*A^3*b*c^7*d^4*e^4 -
 (1215*A^3*b^4*c^4*d*e^7)/4 - 1296*A*B^2*b^2*c^6*d^5*e^3 + 2592*A*B^2*b^3*c^5*d^4*e^4 - 1782*A*B^2*b^4*c^4*d^3
*e^5 + (999*A*B^2*b^5*c^3*d^2*e^6)/2 - 5832*A^2*B*b^2*c^6*d^4*e^4 + 4698*A^2*B*b^3*c^5*d^3*e^5 - (3267*A^2*B*b
^4*c^4*d^2*e^6)/2 - (405*A*B^2*b^6*c^2*d*e^7)/8 + 2592*A^2*B*b*c^7*d^5*e^3 + (1809*A^2*B*b^5*c^3*d*e^7)/8)/b^1
2 - (3*(-c*(b*e - c*d))^(1/2)*(((d + e*x)^(1/2)*(234*A^2*b^4*c^3*e^6 + 4608*A^2*c^7*d^4*e^2 + 9*B^2*b^6*c*e^6
+ 6624*A^2*b^2*c^5*d^2*e^4 + 1152*B^2*b^2*c^5*d^4*e^2 - 1728*B^2*b^3*c^4*d^3*e^3 + 864*B^2*b^4*c^3*d^2*e^4 - 9
0*A*B*b^5*c^2*e^6 - 9216*A^2*b*c^6*d^3*e^3 - 2016*A^2*b^3*c^4*d*e^5 - 144*B^2*b^5*c^2*d*e^5 - 4608*A*B*b*c^6*d
^4*e^2 + 1152*A*B*b^4*c^3*d*e^5 + 8064*A*B*b^2*c^5*d^3*e^3 - 4896*A*B*b^3*c^4*d^2*e^4))/(8*b^8) - (3*(-c*(b*e
- c*d))^(1/2)*((3*A*b^12*c^2*e^5 - 24*A*b^11*c^3*d*e^4 + 9*B*b^12*c^2*d*e^4 + 24*A*b^10*c^4*d^2*e^3 - 12*B*b^1
1*c^3*d^2*e^3)/b^12 - (3*(64*b^11*c^2*e^3 - 128*b^10*c^3*d*e^2)*(-c*(b*e - c*d))^(1/2)*(d + e*x)^(1/2)*(16*A*c
^3*d^2 - B*b^3*e^2 + 5*A*b^2*c*e^2 - 8*B*b*c^2*d^2 - 20*A*b*c^2*d*e + 8*B*b^2*c*d*e))/(64*b^8*(b^5*c^2*d - b^6
*c*e)))*(16*A*c^3*d^2 - B*b^3*e^2 + 5*A*b^2*c*e^2 - 8*B*b*c^2*d^2 - 20*A*b*c^2*d*e + 8*B*b^2*c*d*e))/(8*(b^5*c
^2*d - b^6*c*e)))*(16*A*c^3*d^2 - B*b^3*e^2 + 5*A*b^2*c*e^2 - 8*B*b*c^2*d^2 - 20*A*b*c^2*d*e + 8*B*b^2*c*d*e))
/(8*(b^5*c^2*d - b^6*c*e)) + (3*(-c*(b*e - c*d))^(1/2)*(((d + e*x)^(1/2)*(234*A^2*b^4*c^3*e^6 + 4608*A^2*c^7*d
^4*e^2 + 9*B^2*b^6*c*e^6 + 6624*A^2*b^2*c^5*d^2*e^4 + 1152*B^2*b^2*c^5*d^4*e^2 - 1728*B^2*b^3*c^4*d^3*e^3 + 86
4*B^2*b^4*c^3*d^2*e^4 - 90*A*B*b^5*c^2*e^6 - 9216*A^2*b*c^6*d^3*e^3 - 2016*A^2*b^3*c^4*d*e^5 - 144*B^2*b^5*c^2
*d*e^5 - 4608*A*B*b*c^6*d^4*e^2 + 1152*A*B*b^4*c^3*d*e^5 + 8064*A*B*b^2*c^5*d^3*e^3 - 4896*A*B*b^3*c^4*d^2*e^4
))/(8*b^8) + (3*(-c*(b*e - c*d))^(1/2)*((3*A*b^12*c^2*e^5 - 24*A*b^11*c^3*d*e^4 + 9*B*b^12*c^2*d*e^4 + 24*A*b^
10*c^4*d^2*e^3 - 12*B*b^11*c^3*d^2*e^3)/b^12 + (3*(64*b^11*c^2*e^3 - 128*b^10*c^3*d*e^2)*(-c*(b*e - c*d))^(1/2
)*(d + e*x)^(1/2)*(16*A*c^3*d^2 - B*b^3*e^2 + 5*A*b^2*c*e^2 - 8*B*b*c^2*d^2 - 20*A*b*c^2*d*e + 8*B*b^2*c*d*e))
/(64*b^8*(b^5*c^2*d - b^6*c*e)))*(16*A*c^3*d^2 - B*b^3*e^2 + 5*A*b^2*c*e^2 - 8*B*b*c^2*d^2 - 20*A*b*c^2*d*e +
8*B*b^2*c*d*e))/(8*(b^5*c^2*d - b^6*c*e)))*(16*A*c^3*d^2 - B*b^3*e^2 + 5*A*b^2*c*e^2 - 8*B*b*c^2*d^2 - 20*A*b*
c^2*d*e + 8*B*b^2*c*d*e))/(8*(b^5*c^2*d - b^6*c*e))))*(-c*(b*e - c*d))^(1/2)*(16*A*c^3*d^2 - B*b^3*e^2 + 5*A*b
^2*c*e^2 - 8*B*b*c^2*d^2 - 20*A*b*c^2*d*e + 8*B*b^2*c*d*e)*3i)/(4*(b^5*c^2*d - b^6*c*e)) + (atan((((((d + e*x)
^(1/2)*(234*A^2*b^4*c^3*e^6 + 4608*A^2*c^7*d^4*e^2 + 9*B^2*b^6*c*e^6 + 6624*A^2*b^2*c^5*d^2*e^4 + 1152*B^2*b^2
*c^5*d^4*e^2 - 1728*B^2*b^3*c^4*d^3*e^3 + 864*B^2*b^4*c^3*d^2*e^4 - 90*A*B*b^5*c^2*e^6 - 9216*A^2*b*c^6*d^3*e^
3 - 2016*A^2*b^3*c^4*d*e^5 - 144*B^2*b^5*c^2*d*e^5 - 4608*A*B*b*c^6*d^4*e^2 + 1152*A*B*b^4*c^3*d*e^5 + 8064*A*
B*b^2*c^5*d^3*e^3 - 4896*A*B*b^3*c^4*d^2*e^4))/(8*b^8) - (3*((3*A*b^12*c^2*e^5 - 24*A*b^11*c^3*d*e^4 + 9*B*b^1
2*c^2*d*e^4 + 24*A*b^10*c^4*d^2*e^3 - 12*B*b^11*c^3*d^2*e^3)/b^12 - (3*(64*b^11*c^2*e^3 - 128*b^10*c^3*d*e^2)*
(d + e*x)^(1/2)*(A*b^2*e^2 + 16*A*c^2*d^2 - 8*B*b*c*d^2 + 4*B*b^2*d*e - 12*A*b*c*d*e))/(64*b^13*d^(1/2)))*(A*b
^2*e^2 + 16*A*c^2*d^2 - 8*B*b*c*d^2 + 4*B*b^2*d*e - 12*A*b*c*d*e))/(8*b^5*d^(1/2)))*(A*b^2*e^2 + 16*A*c^2*d^2
- 8*B*b*c*d^2 + 4*B*b^2*d*e - 12*A*b*c*d*e)*3i)/(8*b^5*d^(1/2)) + ((((d + e*x)^(1/2)*(234*A^2*b^4*c^3*e^6 + 46
08*A^2*c^7*d^4*e^2 + 9*B^2*b^6*c*e^6 + 6624*A^2*b^2*c^5*d^2*e^4 + 1152*B^2*b^2*c^5*d^4*e^2 - 1728*B^2*b^3*c^4*
d^3*e^3 + 864*B^2*b^4*c^3*d^2*e^4 - 90*A*B*b^5*c^2*e^6 - 9216*A^2*b*c^6*d^3*e^3 - 2016*A^2*b^3*c^4*d*e^5 - 144
*B^2*b^5*c^2*d*e^5 - 4608*A*B*b*c^6*d^4*e^2 + 1152*A*B*b^4*c^3*d*e^5 + 8064*A*B*b^2*c^5*d^3*e^3 - 4896*A*B*b^3
*c^4*d^2*e^4))/(8*b^8) + (3*((3*A*b^12*c^2*e^5 - 24*A*b^11*c^3*d*e^4 + 9*B*b^12*c^2*d*e^4 + 24*A*b^10*c^4*d^2*
e^3 - 12*B*b^11*c^3*d^2*e^3)/b^12 + (3*(64*b^11*c^2*e^3 - 128*b^10*c^3*d*e^2)*(d + e*x)^(1/2)*(A*b^2*e^2 + 16*
A*c^2*d^2 - 8*B*b*c*d^2 + 4*B*b^2*d*e - 12*A*b*c*d*e))/(64*b^13*d^(1/2)))*(A*b^2*e^2 + 16*A*c^2*d^2 - 8*B*b*c*
d^2 + 4*B*b^2*d*e - 12*A*b*c*d*e))/(8*b^5*d^(1/2)))*(A*b^2*e^2 + 16*A*c^2*d^2 - 8*B*b*c*d^2 + 4*B*b^2*d*e - 12
*A*b*c*d*e)*3i)/(8*b^5*d^(1/2)))/(((135*A^3*b^5*c^3*e^8)/8 - 1728*A^3*c^8*d^5*e^3 - 3996*A^3*b^2*c^6*d^3*e^5 +
 1674*A^3*b^3*c^5*d^2*e^6 + 216*B^3*b^3*c^5*d^5*e^3 - 378*B^3*b^4*c^4*d^4*e^4 + 216*B^3*b^5*c^3*d^3*e^5 - (189
*B^3*b^6*c^2*d^2*e^6)/4 + (27*A*B^2*b^7*c*e^8)/32 + (27*B^3*b^7*c*d*e^7)/8 - (243*A^2*B*b^6*c^2*e^8)/32 + 4320
*A^3*b*c^7*d^4*e^4 - (1215*A^3*b^4*c^4*d*e^7)/4 - 1296*A*B^2*b^2*c^6*d^5*e^3 + 2592*A*B^2*b^3*c^5*d^4*e^4 - 17
82*A*B^2*b^4*c^4*d^3*e^5 + (999*A*B^2*b^5*c^3*d^2*e^6)/2 - 5832*A^2*B*b^2*c^6*d^4*e^4 + 4698*A^2*B*b^3*c^5*d^3
*e^5 - (3267*A^2*B*b^4*c^4*d^2*e^6)/2 - (405*A*B^2*b^6*c^2*d*e^7)/8 + 2592*A^2*B*b*c^7*d^5*e^3 + (1809*A^2*B*b
^5*c^3*d*e^7)/8)/b^12 - (3*(((d + e*x)^(1/2)*(234*A^2*b^4*c^3*e^6 + 4608*A^2*c^7*d^4*e^2 + 9*B^2*b^6*c*e^6 + 6
624*A^2*b^2*c^5*d^2*e^4 + 1152*B^2*b^2*c^5*d^4*e^2 - 1728*B^2*b^3*c^4*d^3*e^3 + 864*B^2*b^4*c^3*d^2*e^4 - 90*A
*B*b^5*c^2*e^6 - 9216*A^2*b*c^6*d^3*e^3 - 2016*A^2*b^3*c^4*d*e^5 - 144*B^2*b^5*c^2*d*e^5 - 4608*A*B*b*c^6*d^4*
e^2 + 1152*A*B*b^4*c^3*d*e^5 + 8064*A*B*b^2*c^5*d^3*e^3 - 4896*A*B*b^3*c^4*d^2*e^4))/(8*b^8) - (3*((3*A*b^12*c
^2*e^5 - 24*A*b^11*c^3*d*e^4 + 9*B*b^12*c^2*d*e^4 + 24*A*b^10*c^4*d^2*e^3 - 12*B*b^11*c^3*d^2*e^3)/b^12 - (3*(
64*b^11*c^2*e^3 - 128*b^10*c^3*d*e^2)*(d + e*x)^(1/2)*(A*b^2*e^2 + 16*A*c^2*d^2 - 8*B*b*c*d^2 + 4*B*b^2*d*e -
12*A*b*c*d*e))/(64*b^13*d^(1/2)))*(A*b^2*e^2 + 16*A*c^2*d^2 - 8*B*b*c*d^2 + 4*B*b^2*d*e - 12*A*b*c*d*e))/(8*b^
5*d^(1/2)))*(A*b^2*e^2 + 16*A*c^2*d^2 - 8*B*b*c*d^2 + 4*B*b^2*d*e - 12*A*b*c*d*e))/(8*b^5*d^(1/2)) + (3*(((d +
 e*x)^(1/2)*(234*A^2*b^4*c^3*e^6 + 4608*A^2*c^7*d^4*e^2 + 9*B^2*b^6*c*e^6 + 6624*A^2*b^2*c^5*d^2*e^4 + 1152*B^
2*b^2*c^5*d^4*e^2 - 1728*B^2*b^3*c^4*d^3*e^3 + 864*B^2*b^4*c^3*d^2*e^4 - 90*A*B*b^5*c^2*e^6 - 9216*A^2*b*c^6*d
^3*e^3 - 2016*A^2*b^3*c^4*d*e^5 - 144*B^2*b^5*c^2*d*e^5 - 4608*A*B*b*c^6*d^4*e^2 + 1152*A*B*b^4*c^3*d*e^5 + 80
64*A*B*b^2*c^5*d^3*e^3 - 4896*A*B*b^3*c^4*d^2*e^4))/(8*b^8) + (3*((3*A*b^12*c^2*e^5 - 24*A*b^11*c^3*d*e^4 + 9*
B*b^12*c^2*d*e^4 + 24*A*b^10*c^4*d^2*e^3 - 12*B*b^11*c^3*d^2*e^3)/b^12 + (3*(64*b^11*c^2*e^3 - 128*b^10*c^3*d*
e^2)*(d + e*x)^(1/2)*(A*b^2*e^2 + 16*A*c^2*d^2 - 8*B*b*c*d^2 + 4*B*b^2*d*e - 12*A*b*c*d*e))/(64*b^13*d^(1/2)))
*(A*b^2*e^2 + 16*A*c^2*d^2 - 8*B*b*c*d^2 + 4*B*b^2*d*e - 12*A*b*c*d*e))/(8*b^5*d^(1/2)))*(A*b^2*e^2 + 16*A*c^2
*d^2 - 8*B*b*c*d^2 + 4*B*b^2*d*e - 12*A*b*c*d*e))/(8*b^5*d^(1/2))))*(A*b^2*e^2 + 16*A*c^2*d^2 - 8*B*b*c*d^2 +
4*B*b^2*d*e - 12*A*b*c*d*e)*3i)/(4*b^5*d^(1/2))